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Abstract 

The paper presents the latest progress on the searches for Lepton Number 
Violating (LNV) B Meson decays, the Lepton Flavour Violating (LFV) decay 
r~ —7- fj,^ fj,~ fi'^ , and the Lepton and Baryon Number Violating (LNV and BNV) 
decays r~ — >• //"'"/i~p and t~ — )• /i~//~p at the LHCb. These searches have 
been performed at a hadron collider for the first time. In the absence of signal 
we put upper limits, which are as follows: B{t^ — ?■ fi^ fj.^ fi^) < 6.3 x 10~*, 
B{t- /i~/u+p) < 3.4 X 10"^ B{t- /i^^^p) < 4.4 x lO"'^ at 90% CL. 
Inclusion of charged conjugate processes are implied throughout this document. 

1 Introduction 

Lepton Flavour Violation has long been observed in the neutrino sector [Ij; the phe- 
nomenon is known as neutrino oscillation. Charged LFV also arises in the Standard 
Model (SM) [2] from the neutrino mass terms via neutrino oscillation at loop level. 
This effect is suppressed by powers of m^/m^, meaning the branching ratios are typi- 
cally < 10~^^|j3j which is well below the observable level. However, many New Physics 
(NP) models (such as MSSM, R-party violating Supersymmetry and Littlest Higgs 
with T-parity |1]) predict LFV decays at much higher rates, significant enough to fall 
within the current experimental sensitivity in certain regions of the model parameter 
space. Whereas any direct observation of charged LFV will be a clear indication of 
NP, exclusion limits for LFV also serve as a powerful tool to exclude parameter spaces 
of NP models. 

Many NP models also feature a wider class of processes: Lepton Number Violation 
(LNV) [5] and Baryon Number Violation (BNV) which are strictly forbidden in the 
SM; most of them [6j predict \B — L\ = 0, 2 where B is the baryon number and L is 
the lepton number. Whereas BNV is highly speculated since it is one of the Sakharov 
[7] conditions to explain the matter-antimatter asymmetry, LNV has not yet been 
observed despite over 70 years of extensive searches in neutrinoless double /3 decay 
[8]. 



2 Majorana Neutrino and LNV in B Meson Sector 

Neutrinoless hadron LNV decays with like sign dilepton final state provide an im- 
portant probe for the existence of Majorana neutrinos. In Majorana neutrino model, 
LNV B meson decays B~ — )■ are produced via two mechanisms at the lowest 

order, which involve either an on-shell or virtual Majorana neutrino (analogous to 
double P decay). On-shell production of — > 7i^fi~fi~ and B~ — )■ D^^u^yti" enables 
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(a) Virtual Majorana neutrino (b) On shell Majorana neutrino 



Figure 1: Example of the lowest order Feyman diagrams for B meson LNV decays 
via Majorana neutrino. 

us to probe Majorana neutrinos in mass range up to 5140 MeV/c^ Beyond this mass 
the modes B^ — )■ D+/i~/i~ and B~ — )■ D+*yU~/i~ are more restrictive. Figure 2 shows 




Figure 2: Lowest order Feyman diagram for the B — )■ tt^D /i decay. 

the four-body decay B~ — )■ 7r~^ fj," jj^ , which was first analysed at the LHCbfQl, fTO]. 
In this case the accessible mass of Majorana neutrinos is smaller, between 260 MeV 
and 3300 GeV, but the rate is enhanced by W coupling. 

In particular, for LNV decays occuring via a fourth massive Majorana neutrino 
z/4, the observation of a LNV decay can provide not only information on the mass 
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m,y^ but also the W,y4i coupling strength \V4i\. In absence of signal in all analysed 
LNV decays, limits on |V4^| coupling were set, which are presented in Figure 3. The 
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Figure 3: Exclusions on Majorana neutrino coupling as function of neutrino mass m^, 
for (a) TT and (b) B~ Dj^/i^/i^ 

current limits for various LNV decays from B factories and LHCb are summarized in 
Table 1. 

3 LFV, LNV and BNV in r sector 
3.1 LFV in r sector 

Studies of LFV in r decays has been performed extensively at the B factories due 
to their high efficiency and clean environment. The most studied channels in the B 
factories are r — t- 3/i and r — > fi'-f; their current experimental limits are 3.3 x 10^® 
from BaBar [17j and 2.1 x lO"*^ from Belle jH], and 4.4 x 10"^ from BaBar jl^ and 
4.5 X 10-8 from Belle |20], respectively (all at 90% CL). 

We focus on the channel t~ — )■ fi'^fi~fi~. In the SM it has a branching ratio smaller 
than 10"^^. Historically, this channel was studied in the B factories, where r are 
produced in pairs in a clean environment {e^e" — ?■ r^r^). Using the thrust axis, effi- 
cient geometry tag of the other r can be performed thereby reducing the combinatorial 
background and providing direct measurement of the number of r produced. At the 
LHCb, the dominant mode for r production (78%) is the leptonic decay D" — )■ t~Vj. 
and r tagging is not possible. This poses big experimental challenges to searches for 
decays like — )■ fi^ fi^fi^ at the LHCb. Nevertheless, the inclusive r cross section 
of 79.5 ± 8.3 yub is large at the LHCb compared to 0.919 nb at the B factories; this 
means two orders of magnitude more of r leptons are produced in the LHCb in one 
nominal year than in the entire run of the B factory experiments. In addition, final 
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Table 1: Current limits on lepton number violating charm (a) and bottom (b) meson 
decays. 
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(b) bottom decays 
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state muons have clean detector signatures thus studies of r jj, ji jjL^ and similar 
decays are totally viable at the LHCb. 

3.2 Search Strategy 

The search is performed by excluding a region of ±30MeV/c2 around the r mass 
until all the analysis choices are finalized (blind analysis). After passing the trigger, 
events are selected using loose cuts based on the kinematics of the reconstructed par- 
ticles. Candidate events are then classified in a three-dimensional likelihood space: 
two multivariate classifiers A^sbody and A^pid, and the invariant mass of the r candi- 
date. A^sbody distinguishes displaced 3-body decays from N(> 3)-body and separate 
combinations of tracks from different vertices using the kinematic and geometrical 
properties of the r candidate. A^pid quantifies the compatibility of each of the three 
decay products with the muon hypothesis using information from the RICH detector, 
calorimeters and muon chambers. Both classifiers are trained on signal and inclusive 
hh and cc background MC, and calibrated with the control channels — )■ (j){fifi)TT~^ 
and J/tp —> fifi for Alsbody and Alpio respectively. The space of each classifier is 
then binned; the number and boundaries of the bins are optimized using CLs[22j 
method. In both cases the optimal number of bins has been found to be 5. For the 
invariant mass classification, the signal mass window of ±15MeV/c^ the expected r 
mass is divided into 6 equally spaced bins and the signal shape is taken from the fit to 
— )■ 0(/iyu)7r+. Both the central value of the mass window and the mass resolution 
are then corrected using the measured scaling and resolution at the LHCb. Figure 4 




(a) Distribution for simulated background (b) Distribution for simulated background 
and the simulated signal as a function of and the simulated signal as a function of 
the PID classifier. the 3 body classifier. 



Figure 4: Distribution of signal events in the two multivariate likelihoods for signal 
(blue / solid) and background (red / dashed) 

shows the distribution of the background and — > fi^fi^fi^ signal MC candidates 
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in each classifier, along with the binning. The vertical lines in Figure 4 are the bin 
boundaries. To calculate the branching ratio, the number of observed signal events 
is normalized to the number of events in the calibration channel — t- 0(yU+/i~)7r~ 



^ _ i3(D+ ^ )7r+) £norm N ,ig 



Norm 

assuming negligible contribution from non-resonant events as suggested by data. 
f{Ds) is the fraction of r~ produced from Dg decays; this factor is required since 
not all r leptons are produced from — t~u^; it is determined using the bb and 
cc cross sections and the inclusive b ^ t~ and c — )■ r~. e^jg and enorm denote the 
total efficiencies for the signal and normalization channels, which takes generation, 
selection and trigger efficiencies into account. 
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(b) Fit to observed events for r~ — > ^+/i^p 
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(c) Fit to observed events for r ^• 
Figure 5: Fit to the events observed in highest likelihood bins for r — t- fififi and 

r — ^ yUyUp 



6 



3.3 LNV and BNV in r Sector 



We focus on the decays — > pfi^ fi'^ and r~ — )■ pfi^fi'^ [23j . Both decays have 
|i? — L| =0 which is predicted by many NP models. The analysis for these channels 
follow closely that of the — t- fi^fi^fi^ mode as described above. The main difference 
here is instead of PID BDT, hard PID cuts are applied on the muon and proton delta 
log likelihoods, which are optimized on signal MC and outer data side bands. Due to 
this hard cut, the normalization factor A^norm is larger. 

4 Results 

The expected number of background events in the signal mass region is determined 
by interpolating from the invariant mass sidebands regions. In case of r — t- fififi 
the peaking SM background Dj" — )■ ri[fifi'-f)fi~^i>f^ is taken into account by the fit, as 
this has been found to be the most relevant exclusive background. Combinatorial 
background is modelled in both cases with an exponential function. The shapes of 
the invariant mass and the multivariate classifier response of the r7(/x/x7)/i^t'^ 
decay are determined with a Monte Carlo sample corresponding to 5fb^^ of data. For 
the r — )■ p/i/i analysis, the expected number of background events is determined by 
interpolating from the mass sidebands also with an exponential and linear function. 
The difference is taken as systematic uncertainty. In the absence of signal an upper 
limit is calculated using CLs method. The results are listed in Table 2. 

Table 2: Limits on the branching fraction obtained by LHCb. 

Decay Limit CL 

T- 12-12-12- <6.3xl0"« @90%CL 

T- i2+f2-p < 3.4 X 10-^ @90% CL [23] 

T-^f2-f2-p < 4.4x10-^ @90%CL |23] 



5 Conclusions 

The LHCb performed its first measurements of LFV r decays. The present limit for 
the r~ — )■ f2~fi~ f2~ channel from the LHCb is still a factor 3 less restrictive than the 
ones set by the B factories; however after the LHCb upgrade the expected limit with 
full data sample is foreseen to be around 8 x 10~^. The LHCb also performed the 
first searches for the LNV and BNV channels: — )■ f2'f2~p, — )■ fi~^f2~p. 
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